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ABSTRACT 

We present the results of a comparative study of the rest-frame optical and rest-frame ultraviolet 
morphological properties of 117 star-forming galaxies (SFGs), including BX, BzK, and Lyman break 
galaxies with B < 24.5, and 15 passive galaxies in the region covered by the Wide Field Camera 3 
Early Release Science program. Using the internal color dispersion (ICD) diagnostic, we find that the 
morphological differences between the rest- frame optical and rest-frame UV light distributions in 1.4 < 
z < 2.9 SFGs are typically small (ICD ~ 0.02). However, the majority are non-zero (56% at > 3cr) 
and larger than we find in passive galaxies at 1.4 < z < 2, for which the weighted mean ICD is 0.013. 
The lack of morphological variation between individual rest-frame ultraviolet bandpasses in z ~ 3.2 
galaxies argues against large ICDs being caused by non-uniform dust distributions. Furthermore, the 
absence of a correlation between ICD and galaxy UV-optical color suggests that the non-zero ICDs in 
SFGs are produced by spatially distinct stellar populations with different ages. The SFGs with the 
largest ICDs (> 0.05) generally have complex morphologies that are both extended and asymmetric, 
suggesting that they are mergers-in-progress or very large galaxies in the act of formation. We also 
find a correlation between half-light radius and internal color dispersion, a fact that is not reflected 
by the difference in half-light radii between bandpasses. In general, we find that it is better to use 
diagnostics like the ICD to measure the morphological properties of the difference image than it is to 
measure the difference in morphological properties between bandpasses. 

Subject headings: cosmology: observations — galaxies: formation - galaxies: high-redshift - galaxies: 
structure 



1. INTRODUCTION 



The Hubble Sequence of 
is seen out to z ~ 1.5 
van den Bergh et"aD 



Brinchmann et al.l 1998 



1991 



Lilly et al 



ralaxies (iHubbld fl936) 
Glaz ebrook et all 1995; 
Griffiths et all UMQt. 
1998; Simard et all 



1999T Ivan Dokkum et all 120001: IStaiiford et al. I built 

Ravindranath et all 12004 ) beyond which high-redshift 
galaxies typically appear compact and irregular (e.g 



isco et allll996t iLowenthal et allll997t iDickinsonl 
van den Berghl 120011 ). Galaxies at high redshift 
were i nitially identified u sing the Lyman-break tech- 
nique (jSteidel et all I1996D . wherein z > 2.5 galaxies 
are selected using a flux discontinuity in the continuum 
(e.g. a U-band dropout) caused by absorption from 
intervenin g neu t ral hydrogen (e.g. iRafelski et all 1 2009 1: 
Yan et all 120091: lOesch etail 12010 : iBunker et all 1200a 
Hathi et alll2010D Subsequent techniques, such as BzK 
(|Daddi et al.l I2004D and 'BX' (jAdelberger et all l200l 
color sele ction, as well as narrow -band selection of Lya 
emitters (|Hu fc McMahonl [1996) , have allowed for the 
photometric identification of star-forming galaxies over 
a wider redshift range. 

The morphological properties of high-redshift galaxies 
are quantified using a wide range of tech niques, including 
the concentration and asymettry indices (|Conselice et all 



120051) . Gini coe fficients (ILotz et a l. 2006), and Sersic pro- 
file fitting (e.g. iCassata et all 120051 ). among others. At 
2 < z < 3.5, star-forming galaxy (SFG) sizes range from 
< 1 kpc to ~ 5 kpc, with the largest often exhibiting mul- 
tiple photometric components (e.g. IBouwens et al.H2004t 
IRavindranath et all 120061; lOesch et all 12009ft . Morpho- 
logical analyses have revealed that most of these systems 
are disturbed and disk-like (i.e., with exponential light 
profiles), with only ~ 30% havin g light profiles consis 
tent with galact i c spheroids (e.g., iFerguson et all 1200 
Lotz et al.lT2006t IRavindranath et al.l I2006t IPettv et"al 



2009). Because high-redshift galaxies are typically faint 
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and clumpy/irregular, it can be difficult to quantify their 
morphology in a m eaningful way. Usi ng a sample of 97 
Lyman-a emitters, iBond et all (|2009T ) found that even 
the half-light radius could not be accurately measured 
without S /N > 30 within a 0"6 a perture. 

In 2003, IPapovich etall (pOOl P03) introduced a dif- 
ferential morphological diagnostic known as the inter- 
nal color dispersion (ICD) that could be used to quan- 
tify the morphological differences between bandpasses 
for astronomical objects. When applied to galaxies, this 
proved to be particularly useful for comparing the rest- 
frame ultraviolet light distribution, which often appears 
clumpy even in low-redshift galaxies, to the rest-frame 
optical light distribution. If the young stars in a galaxy 
are distributed differently from the old stars (as is usu- 
ally the case at low redshift), then that galaxy will 
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have a large ICD. Similarly, inhomogeneous dust dis- 
tributions can lead to spatial variations in the level of 
obscuration of rest-frame UV light, leading to a large 
ICD. Their initial study focused on low-redshift galax- 
ies and found that mid-type spirals exhibited the largest 
ICDs (£ ~ 0.2), irregular galaxies were intermediate 
(£ ~ 0.1), and early-type galaxies had negligible mor- 
phological variation between bands (£ ~ 0). The authors 
also concluded that ICDs could be reliably measured for 
high-redshift galaxies so long as the angular resolution 
(> 0.5 beam kpc^ 1 ^ O'.'l at z = 2) and signal-to-noise 
(S/N > 80) were sufficient. 

Foll owing up on this exploratory study, iPapovich et al.l 
(2005) (P05) computed the ICDs of separate flux-limited 
galaxy samples at z ~ 1 and z ~ 2.3 using a combination 
of Hubble Space Telescope (HST) Wide-Field Planetary 
Camera 2 and Near Infrared Camera and Multi-Object 
Spectrometer imaging. The authors found that the ICDs 
of z ~ 1 galaxies were generally larger and more var- 
ied than those of their higher-redshift counterparts, with 
the former sample including many spiral galaxies and 
the latter being dominated by irregulars and interacting 
systems. They interpret these changes as being due to in- 
creasing stellar population heterogeneity at low redshift, 
with gaseous disks forming around older spheroids. In 
addition, they speculated that the high-redshift galaxies 
with large ICDs were actually mergers-in-progress based 
upon their large apparent asymmetries. 

With the installation of the Wide Field Camera 3 
(WFC3) on HST, we now have a new, improved win- 
dow into the rest-frame optical light of galaxies at 
1 S z S 3. The deep NIR images being taken as 
part of the WFC3 Hubble Ultra Deep Field 2009 (GO 
11563: PI Illingworth) and W FC3 Early Release Science 
(ERS, iWindhorst et al.l I2010T ) programs allow for a de- 
tailed study of high-redshift galaxy morphologies down 
to H\qo ~ 26 mag. In this paper, we present the results 
of a comparative study of the rest-frame optical and rest- 
frame ultraviolet morphological properties of 132 SFGs 
in the WFC3 ERS region of GOODS-S, using the ICD 
as our primary diagnostic. Throughout we will use AB 
magnitudes and assume a concordance cosmology with 
H Q = 71 fan s' 1 M pc" 1 , Q m = 0.27, and Sl A = 0.73 
(jSpergel et al.ll2007t ). With these values, 1" = 8.2 physi- 
cal kpc at z = 2.5. 

2. DATA AND SAMPLE 

Our morphological analyses are performed on star- 
forming and passive galaxies at 1.5 < z < 3.6 using 
HST/ ACS and HST/WFC3 imaging in the GOODS-S 
field. The data products and galaxy sample are described 
below. 

2.1. Advanced Camera for Surveys Imaging 

In the Chandra Dee p Field-South, the sou thern half 
of the GOODS survey (jGiavalisco et alj|200l covers ~ 
160 arcmin 2 of sky and has HST/ ACS observations in the 
F435W, F606W, F775W, and F850LP filters (hereafter 
designated B i35 , V 6 o6, ^775, and z 850 ). For this study, 
we use only the /77s image, which was multidrizzled to a 
pixel scale of 60 mas. The effective exposure time of this 
survey is variable across the GOODS area, but for l'/2 
fixed aperture, a typical /775-band, 5 a detection limit is 
m AB = 27.4. 



TABLE 1 
Galaxy Sample Properties 



Sample 


Number 


Redshift Range 


Reference 


SFGmain 


100 


1.4 < z < 2.9 


1 


SFGhighz 


17 


2.9 < z < 3.5 


1 


PassGal 


15 


1.4 < z < 2.0 


2 



References. - (1) Balestra et al. 2010; (2) 
Cameron et al. 2010 



2.2. Wide Field Camera 3 Imaging 

In order to probe the rest-frame optical light of z ~ 2.3 
galaxies, we use HST/WFC3 F160W images (hereafter 
designated Hi 60 images) taken as part of the WFC3 ERS. 
The Hiqq imaging consisted of 20 orbits over ten visits to 
the GOODS-S field (Program 11563), leading to a total 
of 60 exposures. The exposures were then reduced and 
calibrated (Koekemoer et al. 2010, in prep.), incorpo- 
rating SPARS100 dark frames and correcting for resid- 
ual gain differences between the mosaic quadrants. This 
also included the removal of large-scale scattered light 
residuals and satellite trails. The exposures were t hen 
combined using MultiDrizzle (jKoekemoer et al.l2002f) us- 
ing a pixfrac of 0.8 and a square kernel, to produce fi- 
nal drizzled images with a pixel scale of 60 mas. For 
a l'/2 fixed aperture, a typical i/160-band, 5cr detection 
limit is ttiab =27.2. In order to ensure accurate align- 
ment of the IR imaging with the GOODS ACS imaging, 
the WFC3 exposures were individually aligned to the 
GOODS-S v2.0 z-band catalog and to each other. Even 
for single exposures, the number of sources matched to 
the reference catalog was ~ 500 and the final astromet- 
ric solution appears to be robust to ~ 10 mas or better. 
Further details on this reduction are presented in a forth- 
coming paper (Koekemoer et al. 2010, in prep.) 

2.3. Galaxy Sample 

Our galaxy sample is broken down into three subsam- 
ples: 1.4 < z < 2.9 SFGs (SFGmain), 2.9 < z < 3.5 
SFGs (SFGhighz), and 1.4 < z < 2 passive galaxies 
(PassGal). The passive galaxy subsample includes 15 
objects found in the WF C3 ERS region using a YHVz 
color selection technique (|Cameron etalll20iah . No fur- 
ther cuts were made on this subsample. 

The SFG subsamples are drawn from th e VLT/VIMOS 
spectr oscopi c observations descri bed in IPopesso et al.l 
((20091) and IBalestra et all d20Toh . The VLT/VIMOS 
pointings targeted high-redshift galaxies in GOODS-S 
with B < 24.5, in cluding U-band dro pouts, BzK color- 
selected objects (jDaddi et al.l I2004D. sub-U-d r opout s 
(similar to 'BX' selectio n. lAdelberger et all l2004f ). 
and X-ray sources in the iGiacconi et all (|2002T ) and 
iLehmer et all (|2005l ) catalogs. For this analysis, we ex- 
clude objects with redshift quality flag C (poor qual- 
ity) and objects satisfying the BzK color selection for 
passively-evolving galaxies. In order to exclude unob- 
scured AGN from our subsample of SF Gs, we remove 
any ob ject targeted as an X-ray source by IPopesso et alJ 
(2009), as well as any object that is with in 2 " of an 
X-ray sour ce in the e xpand ed catalogs of iVirani et all 
(2006]) and lLuo et ail (I2008I) . 

Of the remaining SFGs in the WFC3 ERS region, 100 
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TABLE 2 

Galaxy Samples 



Number 


Sample 


Q 


5 


2 


447 


SFGmain 


3:32:01.117 


-27:44:04.830 


2.6355 


455 


SFGmain 


3:32:01.302 


-27:42:43.973 


2.2998 


462 


SFGmain 


3:32:01.509 


-27:44:22.574 


2.5119 


479 


SFGmain 


3:32:02.067 


-27:45:13.878 


2.3131 


482 


SFGhighz 


3:32:02.167 


-27:42:30.428 


3.2630 



This table is only a stub. A manuscript with complete tables 
is available at http://www.nicholasbond.com/Bondl004.pdf 

lie within 1.4 < z < 2.9 and are assigned to SFGmain, 
while a further 17 objects are within 2.9 < z < 3.6 and 
are assigned to SFGhighz. This separation is motivated 
by P03, who found that the internal color dispersion is 
most effective at identifying the signatures of heteroge- 
neous stellar populations when the two filters cover op- 
posite sides of the Balmer and 4000 A breaks. Since we 
are using images taken with the Hiqq and /77s filters on 
HST, the light from young stars will dominate both im- 
ages for galaxies at z > 2.9 (that is, those in SFGhighz) 
and we only expect to see a large ICD in one of these 
galaxies if there is a large column of inhomogeneously 
distributed dust. 

The galaxy subsamples are described in Table [TJ For 
more detail on the selection of these subsamples, see the 
references given in the table. The sky positions and red- 
shifts of individual objects in our parent sample are listed 
in Table [2] ordered by catalog number. In this paper, we 
will refer to SFGs according to their line number in ver- 
sion 2.0 of the VIMOS master catalog (jBalestra et al.1 
2010;. excepting the passive galaxies, which will be de- 
noted P1-P15. Although it is listed in Table[2J we exclude 
object 1379 from the majority of our analysis because it 
was not detected in either the -H160 image or the /77s 
image. 

3. METHODOLOGY 

In the following section, we describle our methodology 
for the measurement of fixed-aperture quantities, includ- 
ing half-light radius and light centroid, as well as the 
internal color dispersion, a differential diagnostic of the 
morphological differences between two bandpasses. Our 
ana lysis is based upon techniques presented in P03, P05, 
and (Bond etHI (pOOl . 

3.1. Fixed- aperture magnitudes and half-light radii 

Because star-forming galaxies are often clumpy, con- 
sist of multiple components, and may be interacting with 
nearby galaxies, there is some ambiguity in choosing the 
appropriate aperture for morphological analysis. We will 
use a simple approach in this paper, measuring all quan- 
tities (with the exception of the internal color dispersion, 
see Section |3~2"|) within a fixed l'/2 aperture. 

Our methodology for measuring object magnitudes, 
centroids, and fixed-aperture half-light radii is described 
in detail in lBond et al.l (|2009D . All of these measurements 
are performed on the original drizzled images, prior to 
the point spread function (PSF) convolution described 
below. 

Using Source Extr action software (SExtractor, 
iBertin fc Arnouti 119961 ). we estimate the maximum 
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Fig. 1. — Internal color dispersion as a function of centroid offset 
between ^775 and Higo-band images for a series of Monte Carlo 
simulations of an SFG in the WFC3 ERS, where each point is 
computed with a random realization of the sky noise in the images. 
The dashed line indicates an iterative quadratic fit to the data 
(discarding outliers). 



contamination by interloping sources within the l'/2 
apertures to be ~ 17% based on the background sky 
density of sources with -H160 < 26. However, because 
the SFGs in our sample were selected using ground- 
based photometry, the contamination rate is probably 
much lower because an interloper within ~ 1" of a 
high-redshift galaxy would have blended with it and 
changed its apparent color. Since these "contaminated" 
high-redshift galaxies will sometimes fall outside of the 
high-redshift galaxy color selection regions, they are 
less likely to be selected than their uncontaminated 
counterparts. 

3.2. Internal color dispersion 

The internal color dispersion attempts to quantify the 
differences in the distribution of an object's light be- 
tween two observed bandpasses without making model- 
dependent assumptions. When comparing between the 
rest-UV and rest-optical, large ICDs can arise if the 
young stars in a galaxy are distributed differently from 
the old stars or the dust is inhomogeneously distributed 
within the galaxy, leading to spatial variations in the 
level of obscuration of rest-frame UV light. At low red- 
shift, P03 found that irregular galaxies have £ ~ 0.1, 
while mid- type spiral galaxies have £ ~ 0.2. High- 
redshift star-forming galaxies are morphologically sim- 
ilar to irregular galaxies, but have had less time to form 
previous generations of stars, so the majority would be 
expected to have £ < 0.1. 



4 



Bond et al. 



For a noise-free image, the ICD is given by, 



N 

E(h,i-al lti -l3) 2 

ah,i2) = 1 ^ , (i) 

i=0 

where the sum is performed over N pixels in a chosen 
aperture, and 1^,% are the fluxes in pixel i, a is the 
ratio of the total fluxes between cutout 2 and cutout 1, 
and (5 is the difference in the background between the two 
cutouts. Since we perform SExtract or sky subtractio n 
on all cutouts prior to analysis (see iBond et alJ 12009). 
we will assume = for all of our measurements. 

If we wish to accurately compute a differential mor- 
phological diagnostic between two observed bandpasses, 
we must ensure that (I) the effective PSF is the same 
in the two images and (2) there is no astrometric offset. 
If the PSF is well approximated by a Gaussian in both 
frames, the former requirement can be satisfied by con- 
volving the lower-resolution image with a Gaussian that 
has a = y/al - a\. Unfortunately, the PSF of the WFC3 
H 1 fin -band imaging is not well approximated by a Gaus- 
sian (jBond et al.H2007l ). so we must convolve each image 
with the PSF of the other. We esti mated the PSF of each 
image using a stack of stars from lAltmann et alJ (2006) 
that fell within the WFC3 ERS region and had temper- 
atures consistent with K and M-type stars, a regime in 
which photometric confusion with galaxies is minimal. 
We further restricted our stack to 29 stars that were iso- 
lated (i.e., the only photometric component within 0'.'6), 
unambiguously stellar (SExtractor stellarity > 0.9), and 
unsaturated. Experiments with stars of a variety of spec- 
tral types reveal that the internal color dispersion is in- 
sensitive to the spectral shape of the PSF across the Hi 60 
filter, with ICDs varying by < 10% of the typical ICD 
uncertainty in SFGmain. 

Noise contributions from the sky background in images 
of high redshift galaxies will affect the terms in Equa- 
tion [TJ In order to estimate the amplitude of this con- 
tribution, we extracted 400 5"x5" cutouts from random 
positions on the Hiqq and /77s images. Cutouts were dis- 
carded if they contained any detected sources. We then 
computed a noise term, 

N N 

£ n (a,N) = Median^ /fj - 2a Median(^/i,i h,i) 

i=0 i=0 

JV 

+a 2 Median^/ 2 ,), 

i=0 

(2) 

where the medians were computed over the set of blank 
cutouts. This was then subtracted from the numera- 
tor and denominator of Equation Q] when computing the 
ICDs of the objects in our sample. 

We calculated the ICDs using 0'/6 to l'/5 circular 
apertures on 5"x5" cutouts in the HST/WFC3 and 
HST/ ACS images, each centered at the flux-weighted 
centroid of all SExtractor d etections within 1"2 o f the 
object position s repo rted in iBalestra et alJ (|2010l) and 
ICameron et alJ (|2010[ ). The sizes of the apertures used 
to compute the ICDs for individual galaxies were selected 
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Fig. 2. — Ratio between the uncertainty in internal color disper- 
sion in our Monte Carlo simulations and the uncertainty as de- 
termined by the analytical approximation of P03 (see Equation [4} 
for the objects in our 2 ~ 2.5 SFG subsample. All points have 
a ™si'ms/ cr i>P05 > 1 as a result of correlated noise in the drizzled 
HST images, but the analytical approximation becomes particu- 
larly bad at high signal-to-noise (i.e. small <r^). 

on an object-by-object basis based upon a visual inspec- 
tion of the cutouts. 

3.3. Monte Carlo Simulations 

For this study, we estimate the uncertainties on the 
internal color dispersion measurements for individual 
galaxies in our sample using object-by-object Monte 
Carlo simulations. In order to do this, we first ex- 
tracted the galaxies in our sample from their Hiqq cutout 
using SExtractor (DETECT_MINAREA = 5 and DE- 
TECT.THRESH = 1.65). Mock H wo cutouts were then 
created by adding the extracted light profile to cutouts 
centered at random positions on the Hiqq image (see Sec- 
tion 13. 2|) . We used the same H 160 -extracted galaxy light 
distributions to generate mock J775 cutouts, normaliz- 
ing the light distribution such that the detected flux was 
equal to that in the real J775 image. This accounts for 
the object's /77s — iii60 color. Finally, we added noise 
realizations from the I775 image. 

The other significant source of uncertainty in our mea- 
surements is the unknown astrometric shift between the 
WFC3 and ACS imaging. In order to simulate this, we 
applied random two-dimensional shifts to the mock Hiqq 
cutouts, with a mean amplitude of 10 mas (our astro- 
metric uncertainties, see Section |2"T2"| . 

Our simulations are designed to approximate the scat- 
ter induced by the sky background on ICD measure- 
ments of an object with an intrinsic £ = and will allow 
us to distinguish objects with intrinsic non-zero ICDs. 
Systematic sources of error, such as astrometric offsets 
and non- uniform sky backgrounds, only produce positive 
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scatter, so the error distribution is non-Gaussian and we 
estimated the bias and positive scatter using the median 
and third quartile, where we define, 



of = 1.496[Q3(0 - Median®]. 



(3) 



We performed a total of 400 simulations for each object. 

Figure Q] shows the dependence of the ICD on the rel- 
ative astrometric shift between 2775 and i?i6o for one of 
our simulations. The dashed line indicates an iterative 
quadratic fit to the data (discarding outliers). Each mea- 
surement includes random realizations of the sky noise, 
leading to the scatter seen about the best-fit curve. For 
the objects in our sample, the median astrometric shift 
at which the resulting systematic error in the ICD ex- 
ceeds the random ICD uncertainty is ~ 0.64 pixels, al- 
though the precise behavior depends upon the shape of 
the light distribution. Point sources, for example, will 
exhibit more sensitivity to astrometric errors than well- 
resolved galaxies. 

In FigureEJ we compare the simulated uncertainties for 
the SFGs in our sample to those that would be derived 
from the analytical approximation of P03, 



N 



s(0 



i=0 

N N ' 

E(J2,i-0) a -E(a a ,i-aB M )a 

i=0 i=0 



(4) 



where B\ i and i?2,i are the fluxes in pixel i in a cutout 
extracted from a blank region of sky. As the expres- 
sion shows, large ICD uncertainties are generally found 
when a source is detected at low signal-to-noise in one 
of the images. At high redshift, the morphological varia- 
tions that we see between wavebands are small, even for 
galaxies with a lot of dynamical activity, so sources with 
individual detections that have S/N > 30 can still have 
uncertainties on their ICD comparable to or larger than 
the median ICD in our sample (£ ~ 0.02). As the authors 
of P03 note, their expression will underestimate the un- 
certainties in the presence of correlated noise. Since HST 
images are typically drizzled to correct for oversampling 
of the PSF, our simulations do in fact yield larger uncer- 
tainties than predicted by P03. P05 accounted for this 
discrepancy in drizzled images by correcting the uncer- 
tainty by a single empirically-determined multiplicative 
factor, C(£). Although C(£) = 6 would be a reasonable 
approximation at > 0.01, it would underestimate the 
uncertainty for some of the brighter objects in our sam- 
ple (tripos ?5 0.01) by factors up to ~ 3. Hereafter, we 
will only use the uncertainties derived from our Monte 
Carlo simulations. 

4. RESULTS 

Below we present the results of a morphological analy- 
sis of the high-redshift star-forming and passive galaxies 
in our sample. All measured quantities, including cen- 
troids, magnitudes, half-light radii, and ICDs, of these 
galaxies are compiled in Table [3l 

4.1. Internal color dispersion 
4.1.1. Star-forming galaxies at 1.4 < z < 2.9 

Some example J775 and Hiqq cutouts for 15 SFGs with 
£ > 0.05 are shown in Figured Of the galaxies shown, all 



0.15 - 



~i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 




_l I I I I I I L 



_l I I I I I I I L 



1.5 2 2.5 3 3.5 
Redshift 



"P" 

X 



b 

\ 

o 

CD 



m 



30 



20 



10 - 







4i j± 

X 

.A X * 
X . XXX A 
A A .» X 
j .% ^X " 

V-A& AKA- 



Li_L 



_L 



_L 



1.5 



2 2.5 
Redshift 



3.5 



Fig. 3. — Internal color dispersion (top) and ICD signal- to- noise 
(bottom) as a function of redshift. We indicate 1.4 < z < 2.9 
SFGs with triangles, 2.9 < z < 3.5 SFGs with open squares, and 
1.4 < 2 < 2 passive galaxies with stars. All error bars are drawn 
from our Monte Carlo simulations and galaxies with <r™ > 0.05 are 
not plotted. In the top panel, the dotted line indicates £ = 0.01 and 
in the bottom panel, the two dotted lines indicate an ICD signal- 
to-noise of 1 and 3. Objects with £(F775W, F160W)/(T™ > 3 show 
statistically significant morphological differences between the ^775 
and i/igo-band images. While the majority of star-forming galaxies 
at 1.4 < z < 2.9 lie above this threshold, very few passive galaxies 
or z > 3 galaxies have £ > at the depth and resolution of these 



are in SFGmain except for object 863, w hich was identi- 
fied as a Lya emitter (LAE) at z = 3.1 in Gronwal l et al.1 
(|2007D . While it is possible that the extended V-band 
emission seen in the cutout originates from the same 
high-redshift source as the Lya emission, the large ICD, 
distinct disk-like morphology, and large size compared to 
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the majority of LAEs at that redshift ()Bond et al.ll2009ft 
suggest that it may be a low-redshift galaxy superim- 
posed on a high-redshift LAE. The remaining galaxies 
exhibit complex morphologies, many with multiple com- 
ponents in both the rest-frame optical and rest-frame ul- 
traviolet. Some of the objects, such as 1204 and 1 598, re- 
semble the "chain galaxies" of lCowie et al.l ()1995[ ) , which 
are believed to be massive galaxies in the act of forma- 
tion, while others (447 and 560) have two-component 
morphologies suggestive of mergers. 

The distribution of internal color dispersions is shown 
as a function of redshift in Figure [3] The 100 galaxies 
in SFGmain exhibit a wide range of ICDs, 56 of which 
are non-zero at > 3cr (see bottom panel). Of the galax- 
ies with small uncertainties in the ICD (cr™ < 0.01), 
only 22% have £ < 0.01, suggesting that the majority 
of color-selected SFGs in this redshift range have under- 
lying populations of old stars that are spatially distinct 
from the young stars. The median ICD in this subsample 
is 0.023. 

There is no apparent correlation between £ and z in 
SFGmain, except for a possible downturn at z > 2.6. 
It may be that the ICD becomes less sensitive to stellar 
population heterogeneity as the H wo filter approaches 
the rest-frame 4000 A break. Alternatively, it may be 
due to evolution - galaxies are expected to evolve towards 
less stellar population het erogeneity and low er dust mass 
with increasing redshift (|Haves et al.l 120 10T) . There is a 
small subset of 21 SFGs with a very large ICD (£ > 0.05), 
which appear to make up a much larger fraction of the 
SFG subsample at z < 2. This is consistent with the 
findings of P05, who found that galaxies with very large 
ICDs went from 7% of their if -band-selected sample at 
z ~ 2.3 to 25% at z ~ 1. 

Figure [4] plots £ as a function of /77s — Hi 60 color. 
We find no evidence for a correlation in SFGmain. Al- 
though this is in agreement with the bulk properties of 
the P05 sample, they found that the galaxies with very 
large ICDs had very red /77s — -Hi6o colors. We see no 
evidence for this effect in SFGmain, but considering that 
P05 only had 7 galaxies with £ > 0.05, the apparent 
differences may be due solely to small-number statistics. 
The lack of a correlation between £ and ^775 — -H160 sup- 
ports the hypothesis that the morphological differences 
between the rest-frame ultraviolet and rest-frame opti- 
cal are due primarily to stellar population heterogeneity 
rather than inhomogeneous dust distributions, as greater 
dust columns would both redden the galaxies and pro- 
duce larger ICDs. A large population of old stars can 
also redden a galaxy's /77s — #160 color, but a galaxy 
with a larger fraction of old stars will only have a larger 
ICD if the two populations are spatially segregated and 
the old stars are not dominating the emission in the ^775 
and -H160 bands. 

4.1.2. Passive galaxies 

The ^775 and Hiqq cutouts for the 15 passive galaxies in 
our PassGal subsample are shown in Figure [SJ Objects 
PI, P3, and P15 are faint in ^775 and have er™ > 0.1, 
so there is little we can say about their differential mor- 
phologies. There may be a diffuse component to the rest- 
frame ultraviolet light of P3 that is not present in the 
rest-frame optical, but the ICD is only greater than zero 



at the ~ 2(7 (£ ~ 0.3) level. The ^775 light distribution 
of P14 is also faint {a^ = 0.07), but the ICD is greater 
than zero at ~ 5a. Its rest-frame optical light distribu- 
tion appears disk-like with a compact central component 
that does not appear in the rest-frame ultraviolet image, 
possibly suggesting the presence of a large column of dust 
in the galaxy center. Object P4 is bright in both band- 
passes, but has two components within l'/2 with very dif- 
ferent -Z775 — Hiqq colors, yielding an anomalously large 
ICD. The top left component has color s consistent with 
a passive galaxy ()Cameron et al.l [20100 and the bottom 
right component may be a low-redshift interloper or a 
nearby star-forming galaxy. 

The remaining 10 objects have single-component, com- 
pact morphologies that appear similar in the two band- 
passes. Of the objects with <r™ < 0.01, 4 of 8 also have 
£ < 0.01. For comparison, only 16 of 74 objects in SFG- 
main with er™ < 0.01 have £ < 0.01, suggesting that the 
stellar populations in the passively-selected galaxies are 
more uniformly distributed than those in SFGs. 

4.1.3. Star-forming galaxies at 2.9 < z < 3.5 

For galaxies with 2.9 < z < 3.5, the ^775 and -Hi6o 
bands both cover light blueward of the 4000 A break 
and the ICD becomes a poor tracer of stellar population 
heterogeneity. If the ICDs of the SFGmain subsample 
are dominated by this heterogeneity (as opposed to dust 
inhomogeneity) , we expect to see a sharp dropoff in the 
ICDs of the SFGhighz subsample. This dropoff does ap- 
pear in the top panel of Figure [3l where 8 of 1 1 galaxies 
with cr™ < 0.01 also have £ < 0.01. Of the entire sub- 
sample of 16 galaxies, only 4 have £ > at > 3<r and 
only 1 at > 5a. 

4.2. Half-light radii and centroid offsets 

High-redshift star-forming galaxies sometimes show 
strong morphological variation between bandpasses, as 
indicated by a large ICD, but differential diagnostics like 
the ICD are only useful if they provide us with infor- 
mation that single-band diagnostics, such as half-light 
radius and centroid offset, cannot provide. In Figure [5j 
we plot the rest-frame optical half-light radius as a func- 
tion of the rest-frame UV half-light radius for the SFGs 
in our sample. In most cases, the half-light radii differ 
by < 20%, with the rest- frame UV light typically being 
more extended for the larger objects. The typical SFG 
is marginally larger (~ 5%) in the i?i6o image even after 
PSF convolution. The points in Figure [5] are also color 
coded according to the magnitude of the ICD for that 
object (blue is smallest, red is largest). Although the 
objects with the largest difference in half-light radius do 
generally have large ICDs, there are many objects that 
exhibit a great deal of morphological variation between 
bandpasses with relatively little difference in half-light 
radius. Figure [5] indicates a clear increase in the me- 
dian ICD with increasing rest-frame optical half-light ra- 
dius. Although convolution with the PSF will artificially 
decrease the ICD for objects near the resolution limit 
(r e < 2 kpc), the trend continues to larger radii. The 
largest objects (> 3 kpc) have a median ICD and disper- 
sion that is ~ 2.5 times that of SFGs with 2 < r e < 3 kpc. 
The correlation between ICD and half-light radius may 
indicate that the non-zero ICDs are being caused by star 
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Fig. 4. — Internal color dispersion as a function of ^775 — -H160 
color. Points and error bars are the same as in Figure[3]and galax- 
ies with a™ > 0.05 are not plotted. There is no evidence of a 
correlation between ICD and rest-UV-to-rest-optical color in the 
2.9 < z < 3.5 SFGs (black points). 

formation from gas that is dynamically young; that is, 
gas that is spatially segregated from a dynamically older 
stellar population and has been recently accreted onto 
the galaxy's host dark matter halo. By contrast, clumpy 
star formation that arises from, for example, disk insta- 
bilities or stellar feedback, would not act to increase a 
galaxy's size in the rest-frame optical. 

Unlike half-light radius, the centroid offset does appear 
to give some indication of the morphological variation 
between bandpasses (see Figure [7]). However, there is a 
great deal of scatter, even out to the largest centroid off- 
sets. Galaxies with the largest ICDs (£ > 0.05) typically 
have centroid offsets of ~ O'.'l, corresponding to 0.83 kpc 
at z = 2.3. 

5. DISCUSSION 

The morphological differences between the rest-frame 
optical and rest-frame UV light distributions in 1.4 < 
z < 2.9 SFGs are typically small (£ ~ 0.02), but are 
statistically significant (more than half are non-zero at 
> 3<r) and are larger than we find in the majority of 
passive galaxies at 1.4 < z < 2. We find that these dif- 
ferences are due to heterogeneous stellar populations, in 
agreement with P05, both because the ICDs show no cor- 
relation with UV-optical color and because a sample of 
higher-redshift galaxies shows much smaller ICDs (typi- 
cally £ < 0.01) when compared between two rest-frame 
UV bandpasses. We cannot rule out a contribution to 
the ICD from a non-uniform dust distribution, but these 
two facts suggest that it is sub-dominant. 

The SFGs with the largest internal color dispersion 
(> 0.05, see Figure [8]) generally have complex morpholo- 
gies that are both extended and asymmetric, suggesting 



{(F775W.F160W) 




r F775W ( kpc ) 



Fig. 5. — Half-light radius in ffi6o as a function of half-light ra- 
dius in 7775 for SFGs at 1.4 < z < 2.9, where points are color-coded 
according to their internal color dispersion. The half-light radii are 
similar between the two bandpasses (~ 20% scatter), with a me- 
dian ratio of r^ lsow /r^" 7 ' T5w = 1.02. Objects sometimes show 
strong morphological variation between bandpasses, as indicated 
by a large ICD, but this is not generically accompanied by a large 
difference in half-light radius. 

that they may be mergers-in-progress. This hypothesis 
is difficult to prove with imaging alone, but the fact that 
many of these galaxies appear clumpy in the rest-frame 
optical suggests that it is not just the young stars that 
have inhomogeneous distributions, but also older gener- 
ations of stars that would have had time to mix in a sin- 
gle, steadily-accreting system. Furthermore, P05 showed 
that single galaxies with half-light radii > 3 kpc are rare 
at z ~ 2, and that ~ 50% of the objects in their sample 
with £ > 0.08 are larger than 3 kpc. If the large ICDs 
are simply due to clumpy star formation, it is occurring 
on a scale much larger than is typical for galaxies at this 
redshift. 

Although many of the passive galaxies in our PassGal 
subsample are centrally compact and have ICDs that are 
consistent with zero, in line with the conventional pic- 
ture of this class of galaxies (Ivan Dokk um et all 120041 : 
I Daddi et al.H2005t [Cimatti et al.1120081: iFrarix et aLlbooa 
IDamianov et alj|2009t ISzomoru et aL[|2010D . at least five 
of the galaxies have non-zero ICDs at > 3cr. This 
suggests that a significant fraction of largely passive 
galaxies at this epoch are still a ctively forming stars , 
possibly in clumpy stellar disks (IStockton et al.l 120081 : 
iMcGrath et~aDl2008t ICassata et alJl2010T >. 

It is important that we continue to explore the wave- 
length dependence of high-redshift galaxy morphologies, 
as different regions of the spectrum probe different physi- 
cal components of the host galaxy. Because high-redshift 
galaxies are often clumpy and irregular, diagnostics such 
as the half-light radius and concentration index are good 
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r e (kpc) 

Fig. 6. — Internal color dispersion as a function of fixed-aperture 
half-light radius in Hisq. Points and error bars are the same as 
in Figure \3\ and galaxies with <r™ > 0.05 are not plotted. For 
2.9 < z < 3.5 SFGs, there is a clear increase in the mean ICD as r e 
increases, even for objects much larger than the WFC3 resolution 
limit (r e ~ 1 kpc), possibly indicating that large ICDs are arising 
from systems that are not yet dynamically settled. 



for obtaining a quantitative measure of the morphology 
in a single band but are usually poor at capturing the 
filter-to- filter variations (e.g., see Figure [5]). As a general 
rule for galaxies at z > 2, it is better to measure the mor- 
phological properties of the difference image than the dif- 
ference in morphological properties between bandpasses. 
The internal color dispersion, which gives the root mean 
squared amplitude of the difference image, is one exam- 
ple of this, but any morphological diagnostic that can 
be applied to a single-band light distribution can also 
be applied to a two-band difference image. An obvious 
next step would be to measure higher-order moments 
of the differential light distribution, but given the rela- 
tively small values of £ we find in high-redshift galaxies, 
one would likely require greater imaging depth in both 
bands to draw meaningful conclusions. 
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TABLE 3 
Photometric Properties 
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<5 a 


z 


#160 


^775 


r H h 


r/ b 




a m d 
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-27:42:30.408 


3.263 
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24.971 ±0.103 


0.29 


0.24 


0.022 


0.011 



a Position of WFC3 iii60-band centroid 

b Half-light radius computed within 1'.'2 fixed aperture before convolution 
c Internal color dispersion between F160W and F775W images 

d Uncertainty on internal color dispersion based on Monte Carlo simulations (see Section 1 3 ■ 3 1 
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Fig. 7. — Internal color dispersion as a function of centroid offset 
between the Hiao and ^775 images for a sample of galaxies in the 
WFC3 ERS. Points and error bars are the same as in Figure [3] and 
galaxies with cr™ > 0.05 are not plotted. Although objects are 
more likely to have a large position offset when their ICD is large, 
the two quantities are not tightly correlated and there are many 
objects with £ > 0.01 and an offset < 0"05. 
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Fig. 8. — Hieo and /77s cutouts of 15 1.4 < z < 2.9 star-forming galaxies with £ > 0.05. The cutouts are 1'.'8 (30 pixels) on a side and 
have been convolved with the PSF of the comparison image (i.e. Hieo with 7775 and vice versa). 
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Fig. 9. — -ffi60 and ^"775 cutouts of 15 passive galaxies. The cutouts are 1"8 (30 pixels) on a side and have been convolved with the PSF 
of the comparison image (i.e. Hieo with J775 and vice versa). 



